Stiff monatomic gold wires with a spinning zigzag geometry 
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Using first principles density functional calculations, gold monatomic wires are found to exhibit a 
zigzag shape which remains under tension, becoming linear just before breaking. At room tempera- 
ture they are found to spin, what explains the extremely long apparent interatomic distances shown 
by electron microscopy. The zigzag structure is stable if the tension is relieved, the wire holding its 
chainlike shape even as a free-standing cluster. This unexpected metallic-wire stiffness stems from 
the transverse quantization in the wire, as shown in a simple free electron model. 
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The manipulation of matter at the atomic scale jjj is 
heralding a technological revolution and opening new re- 
search avenues. A spectacular achievement is the recent 
fabrication |^|| of monatomic chains of gold atoms, the 
ultimate thin wires. Metallic nanowire contacts can be 
created with the scanning tunneling microscope J 3 H , 
with mechanically controllable break junctions p|, or 
even with simple tabletop setups J^]. The relationships 
between conduction, geometric, and mechanical proper- 
ties have been studied by simultaneous measurements of 
conductance and applied force by atomistic | M , Q , 
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]^15| model simulations, 
and by first-principles calculations Jl6|— |l9[] . Until very 
recently, however, only indirect experimental informa- 
tion about the structure of the nanocontacts was avail- 
able. This situation changed dramatically after Ohnishi 
et al Q directly visualized nanometric gold wires by 
transmission electron microscopy (TEM). Surprisingly, in 
a bridge of four atoms connecting two gold tips, which 
was stable for more than two minutes, the atoms were 
spaced by 3.5-4.0 A. Later reports pQ] have even in- 
creased this distance up to ~5 A, a value much larger 
than that in Au 2 (2.5 A) and in bulk gold (2.9 A). 
Gold monatomic chains with a length of four or more 
atoms were independently associated by Yanson et al |J 
to the last conductance plateau during stretching (close 
to one conductance quantum 2e 2 /h). The histogram of 
these plateau lengths showed maxima at regular inter- 
vals, which might be related to the distances between 
gold atoms in the wire. 

In this work we study the structure and stability of gold 
monatomic wires by first-principles density-functional 
calculations pl| . We use Siesta p2|, a code designed to 
treat large systems with local basis sets which has been 
already used to study gold clusters |23| . Tests were per- 
formed for Au2 and bulk gold, using both the local den- 
sity approximation (LDA) p4| and the generalized gra- 
dient approximation (GGA) |B5fl, Core electrons were re- 
placed by scalar-relativistic norm-conserving pseudopo- 
tentials . Valence electrons were described with a ba- 



sis set of double-^ s,p and d numerical pseudo-atomic 
orbitals. Real- and reciprocal-space integration grids 
were increased until a total-energy convergence better 
than 2 meV/atom was achieved. The results are in 
very good agreement with previous calculations, using 
the same functionals, and with the experimental geome- 
tries and vibration frequencies p7f . The GGA improves 
the binding and cohesive energies, but not the geome- 
tries, which are the main focus of this work. In the LDA, 
we obtain, for the gold dimer, a bond length Z=2.51 A, 
a vibration frequency v=190 cm -1 , and a binding en- 
ergy D=3.18 eV. For the bulk fee crystal, the calculated 
nearest-neighbor distance, bulk modulus, and cohesive 
energy are d=2.91 A, B=194 GPa, and £ c =4.55 eV 
respectively. In the GGA, the results are 1=2.57 A, 
v=171 cm" 1 , L>=2.72 eV, d=2.98 A, 5=137 GPa, and 
i? c =3.37 eV. The experimental values are i=2.47 A, 
z/=191 cm" 1 , L>=2.29 eV, d=2.87 A, 5=172 GPa, and 
E c =3.78 eV. 

The wire calculations were performed for infinite 
monatomic chains, using periodic boundary conditions, 
as well as for finite wires of various lengths, either free- 
standing or confined between small pyramidal tips. All 
the calculations were repeated with the LDA and the 
GGA, and both ferromagnetic and antiferromagnetic so- 
lutions were searched. In every case, the geometry was 
relaxed until the maximum forces were smaller than 
10 meV/A (16 pN). As an additional cross-check, some 
critical geometries were recalculated with a different 
code, using a plane wave basis set. The results will be 
presented in full elsewhere. In short, we have found no 
qualitative differences, and only very minor quantitative 
differences between the finite and infinite wires, between 
plane wave and local basis sets, and between LDA and 
GGA, and no magnetic solutions could be stabilized at 
any wire length. We present in what follows the Siesta 
LDA results for the infinite wires, except where stated. 

Fig. ^ shows the wire geometry and the binding en- 
ergy as a function of the wire length. Except when very 
stretched, the wire adopts a nonlinear, planar zigzag ge- 
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FIG. 1. First-principles, density-functional results for the 
bond angle a (a), and bond length r (b) in a monatomic gold 
wire with zigzag geometry, as a function of its length per 
atom, (c) Binding energy E in the zigzag (solid symbols) and 
linear wires (open symbols). 

ometry, with two atoms per unit cell. Unconstrained 
relaxations with larger cells did not result in longer peri- 
ods, nor in out-of-plane deformations. The energy shows 
a shallow minimum at a length of 2.32 A/atom, with 
a bond angle of 131°. The stability of this geometry 
was demonstrated by checking that the dynamical ma- 
trix, calculated in a cell of 16 atoms, had no negative 
eigenvalues. For comparison Fig. lc shows the energy of 
a wire constrained to a linear geometry, which has a min- 
imum 0.24 eV/atom higher, and at a wire length 0.25 A 
longer, than in the zigzag geometry. This difference in 
wire length is almost entirely due to the change in bond 
angle, since the bond distances differ by only 0.02 A be- 
tween the two minima. The bond angle increases with 
stretching, but the wire becomes linear only shortly be- 
fore breaking. 

The comparison between the band structures of the 
linear and zigzag wires (Fig. ^) offers some hints for un- 
derstanding their relative stability. In the linear chain, 
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FIG. 2. Electronic band structure of the linear (a) and 
zigzag (b) wires for a length of 2.32 A/atom. The linear-wire 
bands have been folded onto a two-atom Brillouin zone to 
facilitate the comparison. The energies are relative to the 
Fermi level. 

the overlap between the filled d states broadens the d 
bands until they reach the Fermi level, destabilizing the 
wire with their associated high density of states. For the 
same wire length, the zigzag configuration allows a larger 
bond distance, that brings back the d bands below the 
Fermi level and leaves a single s band crossing it. This 
is consistent with the observation of a single conduction 
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FIG. 3. Energy versus wire length for a simple model of 
the wire, considered as a free-electron tube of fixed volume, as 
shown in the inset, with ao = Vab = 3 A and 5 = 2 A. Dashed 
and solid lines correspond to allowing one or two occupied 
bands respectively. 
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FIG. 4. Calculated valence electron density of a 

seven-atom wire suspended between two small pyramidal tips, 
averaged over the three equivalent relaxed configurations, 
with the zigzag plane rotated by 0, 120, and 240° around 
the wire axis. The plot is an isosurface of constant averaged 
pseudo electron density p=0.16 e/bohr 3 . The numbers are 
relaxed distances in angstroms. 



channel in the monatomic wires Q3J. A Peierls dimeriza- 
tion instability is expected since the Fermi wave- vector is 
at the edge of the two-atom Brillouin zone. We have ob- 
served, however, that the magnitude of this gap-opening 
instability is negligible, only slightly noticeable just be- 
fore the wire breaks, and thus playing no substantial role 
in the physics described here. 

Although the appearance of a zigzag instability un- 
der compression may seem natural, its presence in a 
stretched wire is more surprising. Furthermore, its sta- 
bilization at a finite wire length is even harder to under- 
stand, since one would expect the wire to collapse into 
a compact, high-coordination structure typical of met- 
als. However, we find that even free-standing clusters of 
four or eight atoms (the sizes calculated) are also stable 
with a zigzag chain structure. Although unexpected, this 
stability arises very naturally from the transverse quan- 
tization of the electron states. To see this, we model the 
wire as a tube of length a per atom, with a rectangu- 
lar section b x c. Consistently with the standard jellium 
model Jl5| ], we assume a fixed volume per atom abc, but 
we allow a larger 'box' section (b + 8) x (c + 8) to ac- 
count for an electron 'spillage' 8/2 out of each jellium 
edge p8fl . Accepting from the ab-initio calculation that 
the zigzag is planar, we also fix its thickness c or, equiv- 
alently ao — yfab. The resulting free-electron energy is 
shown in Fig. ^ as a function of the wire length a, for rea- 
sonable values of ao and 8. With a single occupied band, 
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FIG. 5. (a) Wire-length dependence of the optical 
phonon frequencies of a zigzag wire at V, calculated using 
the frozen-phonon method. Solid and open circles stand for 
the longitudinal and transversal modes respectively. For com- 
parison, we also show the longitudinal phonon frequency at X 
(which folds to V in a two-atom cell) for a wire constrained to 
be linear (squares), (b) Phonon dispersion curves for a zigzag 
wire of 2.62 A/atom length. 



the compromise between the transversal and longitudi- 
nal kinetic energies results in a single minimum (dashed 
line). Including the second band, which becomes par- 
tially occupied at somewhat shorter lengths, allows the 
energy to decrease again (solid line), reproducing very 
well all the qualitative features observed in the ab-initio 
curve, such as the positions of the maximum, the mini- 
mum, and the point at which the second band crosses the 
Fermi level (1.83 A/atom). The basic physics that this 
model illustrates is the higher stability of certain wire 
sections, due to the transverse quantization of the delo- 
calized electron states [Q. This shell structure effect, 
which has been recently observed for sodium wires p9[ , 
is similar to the so-called magic numbers (particularly 
stable sizes) of small metal clusters pO| , pl[ . The zigzag 
shape is a particular realization of these stable sections 
for the monatomic gold wires. 

In agreement with previous ab initio calculations |l9| l 
we find that the wire becomes unstable and breaks spon- 
taneously when pulled by a force of more than 2.2 nN, 
i.e. beyond a length of 2.9 A/atom, much shorter than 
that apparently observed in stable wires We offer 
here an explanation for this puzzling discrepancy, based 
on the predicted zigzag geometry: if the actual wires ob- 
served have an odd number of atoms, with those at the 
extremes fixed by the contacts, the odd-numbered atoms 
would stay almost fixed on the same axis, while the even- 
numbered ones could rotate rapidly around that axis, of- 
fering a fuzzy image that could be missed by the TEM. 
We have calculated the relaxed geometry and the rotation 
energy barrier for a seven-atom wire suspended between 
two pyramidal tips. We find that the stable geometry 
is almost equal to that of the infinite wire, and that the 
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rotation barrier is only 60 meV for the entire wire. The 
effect is illustrated in Fig. £|, where we show the electron 
density averaged over rotated configurations. Although 
not directly comparable to a TEM image, it can indeed be 
qualitatively appreciated that the odd-numbered atoms 
appear much sharper than the even-numbered ones, giv- 
ing the impression of a four-atom wire with a large inter- 
atomic separation, similar to that observed experimen- 
tally. From the energy barrier obtained, we estimate that 
the thermal rotation would slow down to the millisecond 
scale, allowing the zigzag visualization, only for temper- 
atures below ^40 K. 

Fig. ||(a) shows the calculated transversal and longi- 
tudinal phonon frequencies at T, for the zigzag wire, 
as a function of its length. Negative values indicate 
modes with imaginary frequency, implying the breaking 
of the unstable wire. At the wire's equilibrium length 
(2.32 A/atom), the T-point frequencies are 113 and 
219 cm -1 , for the transversal and longitudinal modes, re- 
spectively. These are quite larger than the bulk phonon 
frequencies, but comparable to those of the dimer. This 
is not surprising if we consider that the wire interatomic 
distance is only slightly larger than that in Au2 . Fig. ||(b) 
shows the phonon dispersion relations for a wire length of 
2.62 A/atom, obtained from the full dynamical matrix in 
a supercell of sixteen atoms, calculated with finite differ- 
ences. H^J We hope that the comparison of the results in 
Fig. ||(a) and (b) with those of point contact spectroscopy 
experiments [ |33| will help to confirm our predicted zigzag 
distorsion. 
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